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Symmetry energy coefficients of asymmetric nuclear matter are investigated as the inverse of nu- 
clear matter polarizabilities with two different approaches. Firstly a general calculation shows they 
may depend on the neutron-proton asymmetry itself. The choice of particular prescriptions for the 
density fluctuations lead to certain isospin (n-p asymmetry) dependences of the polarizabilities. Sec- 
ondly, with Skyrme type interactions, the static limit of the dynamical polarizability is investigated 
C*~) ' corresponding to the inverse symmetry energy coefficient which assumes different values at different 

asymmetries (and densities and temperatures). The symmetry energy coefficient (in the isovector 
channel) is found to increase as n-p asymmetries increase. The spin symmetry energy coefficient is 
also briefly investigated. 
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I. INTRODUCTION 

. The (n-p) symmetry energy coefficient and its dependence on the nuclear density has been extensively studied and 
this is of relevance, for example, for the description of macroscopic nuclear properties as well as for proto-neutron 
and neutron stars. It represents the tendency of nuclear forces to have greater binding energies (E/A) for symmetric 
systems - equal number of protons and neutrons. It contributes as a coefficient for the squared neutron-proton 
asymmetry in usual macroscopic mass formula, E/A = Hq(A, Z) + a T (N — Z) 2 /A 2 , where Hq does not depend on 
"■-^ ■ the asymmetry, Z, N and A are the proton, neutron and mass numbers respectively. Other powers of the asymmetry 
»^ \ (proportional to (N — Z) n for n ^ 2 pj) are usuallyexpected to be less relevant for the equation of state (EOS) of 
nuclear matter based on such parameterizations 0, 13] • The same kind of parameterization is considered for nuclear 
matter where instead of nucleon numbers one has to deal with densities. a T is also the parameter which measures 
the response of the system to a perturbation which tends to separate protons from neutrons. It is given by the 
static polarizability of the system which also may depend on the asymmetry of the medium. This point has been 
developped and emphasized recently 0, 0. The spin symmetry energy coefficient of nuclear matter may also be 
defined, a a , representing the cost in energy to make the system spin-asymmetric (and eventually polarized nuclear 
matter). The spin channel is relevant for the study of the neutrino interaction with matter because it couples with 
axial vector current together with the scalar channel |E0@- A suppression of the spin susceptibility (in this work we 
will be dealing rather with its inverse) leads to the suppression of Gamow Teller transitions which are of interest for 
the supernovae mechanism Q and eventually to instabilities associated to ferromagnetic polarized states 0, . In this 
work we articulate and extend the ideas developped previously for the dependence of symmetry energy coefficients on 
neutron-proton asymmetry. For this we use a calculation for the static polarizabilities - proportional to the inverse of 
the symmetry coefficients in asymmetric matter - which was done using Skyrme effective forces in 0, . 
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II. GENERALIZED SYMMETRY ENERGY COEFFICIENTS 



When considering a small amplitude external perturbation e, the medium polarizability is defined as the ratio of 
the density fluctuation (/3 = Sp n — 5p p ) to the amplitude of the external perturbation and it can be written as |5Hl0|: 

IF-* = I = t (1) 

where A. s ,t is the neutron-proton (isovector) symmetry coefficient (s = 0, t = 1 - spin, isospin) and b = p n /p p — 1 
is an asymmetry coefficient. For the other channels (different (s,t)) one may define different symmetry coefficients. 
Note that A Sl t is a function of b and (3 and these parameters may be related, as argued below. The occurrence of 



these functional dependences of A s ,t can be found just by the first stability condition with respect to the (density) 
fluctuation from which one defines the polarizability QJ: dH/d(3 — 0. 



A. Isospin dependence of As.t 



We consider A s .t to be a function of the density fluctuation (3. Although (3 is not the explicit n-p asymmetry (given 
by 6) we will consider that it depends on it (as it was also argued in |5(). We consider these parameters are related 
to each other and therefore we will write A = A{(3) shortly. In |4( two different prescriptions were discussed for (3 in 
the calculation of the response function of asymmetric nuclear matter. We have used (and it was shown to be the 
more reasonable prescription) the one which leads to the following relation between the fluctuation (3 and the explicit 
asymmetry (6): 

P = 5p n (0), (2) 

Where 5p n is the neutron density fluctuation. In the n-p symmetric limit f3 = 2Sp n and in another limit, in neutron 
matter, (3 = 5p n . The above prescription (expression (J2J) is based on the assumption that the density fluctuations are 
proportional to the respective density of neutrons and protons, i.e., 8p n / (3 = p n /p, being p the total density. In spite 
of being rather well suited for the isovector channel, this kind of assumption can be considered as a starting point 
for the other channels (spin, scalar) in asymmetric nuclear matter. Prescription J5J is therefore model-dependent and 
different choices for it yield other forms for the the (asymmetric) static screening functions. The dynamic response 
functions are less sensitive to this prescription 

From the solution of the polarizability Q we calculate the first derivative with relation to b: 

dp _ ep dA ' 1 _ (3 dA - 1 

db ~ 2A 2 QA db ~ ~AT! db ' (3) 
Another expression can be obtained from the relation between b and /? of (2J ■ It yields: 

d(3 (3 



db (2 + &)(l + &)' 

Equating these two last equations we obtain: 



(4) 



A (2 + 6)(l + 6) " P db ' [ ' 

From which it is possible to derive the following relation between the isospin s.e.c. and the n-p asymmetry [{|: 



A 0,l _ ,0,1 2 + 26 (Q) 



In this expression A sym = a T ~ 30MeV is the s.e.c. of symmetric nuclear matter (6 = 0). For b — 2 (neutron density 
three times larger than the proton density) we obtain A = 1.5A sym . In the limit of neutron matter A(b — > oo) = 

2Asym • 

Other prescriptions for the density fluctuation, leading to different dependences on the n-p asymmetry (6), can be 
given by: 



1) (3 = const — > »4o,i — const, 

2 + b m 2 + 2b m (7) 

2) f3 = 5p n ^ - bm -> Ao,i = Asym 2 - h „ (m > 1 integer). 



The first of these alternative prescriptions leads to a constant symmetry energy coefficient A = a T . However, it 
corresponds to Sp n = Sp p , independently of the asymmetry of the medium. This does not seems to be reasonable for 
example because in the limit of neutron matter there would be no proton density. Furthermore this prescription has 
been used in the dynamical response function yielding seemingly non physical results Q . The second of prescriptions 
@ is a very general one valid for arbitrary integer numbers m and was considered for the sake of simple algebraic 
calculation. Probably m should be not large because it would lead to a too much strong (stiff) dependence on b. 

Another assumption for deriving expressions @ and (JJJ was that p is independent of b. This would be non trivial 
if one considers a complete self consistent calculation with the equation of state of a proto-neutron star, for example. 
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III. POLARIZABILITIES WITH SKYRME FORCES 



A nearly exact expression for the dynamical polarizability of a non relativistic hot asymmetric nuclear matter at 
variable densities was derived with Skyrme interactions in [4| . 

The general static screening function A s j in asymmetric nuclear matter at finite temperature was explicitely written 
m 0. (The coefficient b is related to a frequently used asymmetry coefficient: a — (2pon — Po)/po, by the expression: 
b = 2a/(l - a).) 



A. Results for Skyrme interactions 

The Skyrme interactions used are: SLyb H3, SkSC4, SkSC6 and SkSCIO [HEl. In Figure 1 we show the inverse 
of the static polarizability (generalized symmetry energy coefficient) in the isovector channel for diverse Skyrme forces, 
as a function of the asymmetry parameter b. For symmetric nuclear matter 6 = the different Skyrme forces yield 
values between 27 and 34 MeV, the usual values adopted in the litterature. For increasing b the coefficient increases, 
being the slope strongly dependent on the interaction. We plot one case for dense nuclear matter p = 2po (long 
dashed-short dashed line) with force SLy(b). We note that the increase of Ao,i is much smaller. These results were 
discussed more extensively in Q. 

In Figure 2 the spin symmetry energy is shown as a function of b with forces SLy, SkSC4, SkSC6 and SkSCIO. The 
values for symmetric nuclear matter are very different. The common trend is the increase of with b, i.e., at very 
asymmetric matter the spin interaction tends to become more repulsive. However the particular behavior of the spin 
s.e.c. with b is different for each effective force at a given density. We can compare our results to the ratio of spin 
susceptibility of interacting neutron matter to the non interacting Fermi gas obtained by Fantoni, Sarsa and Schmidt 
|l4j by means of the auxiliary field diffusion Monte Carlo method. This ratio is proportional to the polarizability 
as obtained in expression and therefore inversely proportional to the spin symmetry coefficient ^4i,o- First of all 
we note that, in most cases, the values they find are all positive for the range of densities considered by them, from 
0.75po up to 2.5po (in our calculation the total density was kept constant and equal to the saturation density po). 
The slope seems to be nearly the same as that we obtain for low values of the n-p asymmetry. Consequently they may 
obtain instabilities for higher density neutron matter whereas we do not observe this result in our calculations with 
Skyrme forces at the saturation density (the comparison is meaningful for neutron matter: b very large). A further 
comparison at different total nuclear densities is to be shown elsewhere 



IV. SUMMARY AND CONCLUSIONS 

Summarizing, the dependence of the s.e.c. on the n-p asymmetry was studied extending the results of ref. 
Different prescriptions for the density fluctuations (expressions (\2l7l ) in asymmetric nuclear matter lead to different 
dependences of the symmetry energy coefficients on the asymmetry parameter b. The n-p asymmetry dependence of 
the s.e.c. in the different channels was analyzed for different Skyrme forces. They may yield very different behaviors 
including the possibility (or not) of nuclear matter to undergo phase transitions. These forces can describe different 
behaviors of the symmetry energy coefficients. Therefore, in principle, different values can be expected for the (bulk) 
symmetry energy coefficients in asymmetric nuclear matter with different n-p asymmetries. It would be interesting 
to apply the results obtained here and in [5j in studies of the equation of state of asymmetric nuclear matter, trying 
to extract experimental constraints for Aj\u, in particular in the isospin channel (s = 0, t = 1). (For works which 
can be related with the present ideas see, for example, these proceedings and [l5j). Although we have dealt with 
neutron-proton density asymmetry we can expect that the same kind of ideas can be applied for the neutron-proton 
number asymmetry used in mass formulas for finite nuclei |3j . In the spin channel it is possible to expect spin polarized 
asymmetric matter yielding magnetic fields in neutron stars, as discussed in 0. However with the increase of b we find 
that the spin interaction may be rather repulsive, hindering this magnetization effect with the use of these Skyrme 
forces for the value of density analyzed here (po)- In this work we do not include the simultaneous study of the density 
dependence of the symmetry coefficients. 
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Figure captions 

Figure 1 Neutron-proton symmetry energy coefficient Ao,i = p/(2Il a ^ 1 ) as a function of the asymmetry coefficient 
b for several Skyrme forces: dotted line (SGII), thick dotted (SkSC4), long-short dashed (SLy(b)), crosses (SkSC6), 
thick long-short dashed (SLy(b) at twice p n ), medium thick dotted (SkSCIO). 

Figure 2 Spin symmetry energy coefficient = p/(2Tl)f) as a function of the asymmetry coefficient b for the 
different forces: thick dotted (SkSC4), long-short dashed (SLy(b)), medium thick dotted (SkSC6), thick solid line 
(SkSCIO). 
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